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We perform a systematic study of the thermodynamics of 
quantum gases in the unitarity limit. Our study makes use 
of a "Universality Hypothesis" for the relevant energy scales 
of a many-body system at unitarity. This Hypothesis is sup- 
ported by recent experiments and can be proven in Boltzmann 
regime. It implies a universal form for the grand potential 
which is characterized by only a few universal numbers in 
degenerate limit. This universal form provides a simple way 
to determine the density profile of a trapped fermion super- 
fluid as well as the second sound velocity of a homogeneous 
superfluid at unitarity. 

Feshbach resonance has introduced a whole new di- 
mension in the research of degenerate quantum gases. 
Through this resonance, effective interactions between 
atoms are dramatically increased. Such resonance arises 
when the energy of a pair of scattering atoms is tuned 
close to that of a molecular bound state by an external 
magnetic field, thus causing substantial resonance scat- 
tering. At resonance, scattering reaches the unitarity 
limit with a diverging scattering length a 8C , and a cross 
section reaching the maxium value iir/k 2 , where k is the 
relative momentum of the scattering atoms. These prop- 
erties are universal because they are independent of any 
feature of atomic potentials. 

This universality, simple as it is, poses a challenging 
many-body problem, as there are no small perturbative 
parameters. On the other hand, it can lead to great 
simplification if one makes a very reasonable assumption 
(referred to as "Universality Hypothesis") : that the only 
dominant length scale at unitarity in the ground state is 
inter-particle spacing n -1 / 3 , where n is the density. This 
is based on the idea that the only relevant length scales 
in the ground state are a sc and n -1 / 3 . Since a sc diverges 
at resonance, it must disappear from all physical quan- 
tities, leaving n -1 / 3 the only relevant length scale. The 
word "Hypothesis" is to indicate that although universal- 
ity has emerged in approximate calculations [1] [2], it has 
not been proven rigorously except in Boltzmann regime 
[3] . This hypothesis also implies that for both bosons and 
fcrmions, the only relevant energy scale is the "Fermi" en- 
ergy 8p(n) = (h 2 /2M)(3ir 2 n) 2 / 3 . For the same reason, 
the transition temperature T c of a fermi superfluid must 
scale as Tp = Sp/ks, i.e. Tc = 7?f, where 7 is a uni- 
versal constant. That 7 can be of order 1 was suggested 
by Holland et.al [1]. Current estimates of T c range from 
0.5 to 0.2T F [1] [4]. The possibility of such a high T c 
has made Feshbach resonance a focus of attention in the 
current race for achieving fermion superfluidity. On the 
other hand, the normal state is no less intriguing, for it 
contains the same non-perturbative effects. 

The universal properties of a Fermi gas in unitarity 



regime have been demonstrated recently by a sequence 
of beautiful experiments [5] [6] [7] [8]. John Thomas's 
group [5] has pointed out that that the interaction en- 
ergy Ei n t of a Fermi gas of 6 Li near Feshbach resonance 
is the form j3£ F , f3 « —0.25 over the temperature range 
1 > T/Tp > 0.1. A direct measurement of this energy 
was performed later by Salomon's group who found a sim- 
ilar energy at higher temperatures which remains smooth 
across the resonance. Not only does Ei nt saturate at reso- 
nance despite the diverging scattering length, it scales di- 
rectly with the Fermi energy. Similar saturation are also 
observed in RF spectroscopy by Ketterle's group [7] and 
Jin's group [8] . The observed relation between Ei nt and 
£ F is a support for the Universality Hypothesis. Further 
experiments on other alkali fermions will help to verify 
its validity. The sign of Ei nt and its continuity across the 
resonance [6], however, require additional physics and are 
related to formation of molecules [3] . Since the Duke ex- 
periments cover the temperature range above and below 
the estimated T c , it raises the question of how superflu- 
idity is affected by unitarity, and their signature in the 
unitarity regime. 

In this paper, we perform a systematic study of the 
thermodynamics of quantum gases at unitarity using the 
Universality Hypothesis. We shall show that: (A) At 
unitarity, the thermodynamic potentials acquire univer- 
sal forms which depend only on the nature of the ther- 
modynamic phase. (B) The properties of a degenerate 
Fermi gas near resonance (be it normal or superfluid) 
are characterized by only a few universal numbers. (C) 
Universal thermodynamics provides a simple way to de- 
termine the density profile of a trapped fermion super- 
fluid near resonance. (D) It also allows one to calculate 
the hydrodynamic modes of a fermion superfluid at fi- 
nite temperature. (E) Bose gas in the unitarity limit, if 
stable, will have a fermionic energy density. 

Before proceeding, we first make clear what quantities 
the Universality Hypothesis describes. Let us consider 
the Hemholtz free energy density / = f(T, n, B, {r,}), 
which is a function of T, n, external magnetic field (which 
controls a sc ), and other interaction lengths such as effec- 
tive range, etc, (collectively denoted as {r<}). While a sc 
diverges at resonance (say, at B = B a ), it is assumed that 
{ri} remain atomic size so that Xi — r^n 1 / 3 << 1. If / 
is smooth across the resonance (as indicated in ref. [6]), 
then in the neighborhood of B a it is well approximated 
by its value at resonance f(T,fi;B ,{xi}). Moreover, if 
/ has an asymptotic expansion in x%, we then have 

f(T, n; B , {x,}) = f(T, n; B , 0) (1 + 0( Xi )) . (1) 

Universality Hypothesis describes the first term eq.(l) 
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that it has only two energy scales, £p{n) and fc^T; 

and is independent of B Q for it would otherwise introduce 
additional energy scales. In this paper, quantities of ideal 
gas will be denoted by a superscript "°" . 

I. Universal Thermodynamics: We begin by deriv- 
ing some general thermodynamic properties at resonance 
valid for different phases. To be efficiency, we consider 
the grand potential instead of the Hemholtz free energy. 
The former is defined as fl(T, /zj, V)=— k B T\nTr 
e -(H- MT Ar T - Mi Ar i )/fe B T for a twQ componcnt quantum 

gas (labelled as f and J.) with identical mass M in a 
volume V. We first consider /ij = fi^ which implies 
n-f = ni = n/2, or to = (n-f — njJ/2 = 0. Accord- 
ing to Universality Hypothesis, all microscopic scales are 
absent at resonance. The only energy scales are k B T 
and fi. The corresponding density scales are A~ 3 and 
nf(fj), where A = h/^/2nMk B T is the thermal wave- 
length, and nf(n) = (37r 2 ) -1 (2M/i/?i) 3 / 2 . Since pressure 
P = —Sl/V, we have the following two equivalent rela- 
tions from dimensional analysis 
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- 1 ) = ^W), 
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where x — k B T/u, and (Wo, Go) arc dimensionless scal- 
ing function. These two forms are useful in Boltzmann 
and degenerate regime respectively since their arguments 
are small in these cases. Using the well known relation 
e = Ts + fin — P for energy density e and the Gibbs- 
Duham relation dP = nd/i + sdT, we have 



n f (n) ( G (x) - -xG' (x) 



s = 



2nf{p)ki 



(3) 



(4) 



Eq.(2) to (4) are identical to those of idea gas, where 
k B T and \i are also the only energy scales. However, 
unlike ideal gas, the absence of microscopic energy scale 
here is not due to absence of interaction, but instead 
the maximum scattering between particles. Since uni- 
versality hypothesis makes no reference to the thermo- 
dynamic phase, eq.(2) to (4) apply to both normal and 
superfluid phases, which have of course, different scaling 
functions. The scaling functions, however, are constraint 
by the positivity of s and n, as well as stability conditions 
d 2 P/0T 2 = ds/dT > and d 2 P/diJ 2 = dn/d/j, > 0. 

The density profile n(x) in a non-uniform potential 
V(x) can be readily determined from eq.(3) within local 
density approximation (LDA) by replacing /i with /x(x) = 
}i— V(x). The total energy of a system of N particles E = 
E(T, N) can the be obtained by eliminating fi from the 
relations N = J n(x) = N = N(T, fi) and E = J e(x) = 
/3P(x)/2 = £(T,/i). 

When m 0, Universality Hypothesis implies that 



P(T,u h n L ) 



2un f (u) 



G(x; v/n) 



where /i = (/ij + /UjJ/2, v = fj,^ — fj,^, and Q(x;v/a) is 
a scaling function even in v due to the invariance of f2 
under spin exchange. For small v, we have 

g( X ;u/n) = g (x) + ^-(^j +... (6) 

where Q 2 is another dimensionless function. Defining 
magnetic susceptibility \ and specific heat c at constant /1 
as to — (n-f — ni)/2 = dPjdv = \ v and c = T(ds/dT)^, 
eq.(5) and (6) then implies c/(T X ) = k%Q%{x)/Q 2 {x), 
which is a universal function of x = k B T/n. 

II. Boltzmann limit: This is the limit where the 
fugacities Zi = e AJi//fcisT , i =|,| arc small; and where 
universality can be proved rigorously [3]. We shall first 
generalize the derivation in ref. [3] to arbitrary spin polar- 
ization and the derive the key thermodynamic properties 
in the Boltzmann regime. Expanding £1 in [9] for a 
Fermi gas, we have 

P(T, MT , Mi ) = p(°)(T, MT , Mi ) + 2>/2&2^s T *i (7) 

where P^(T, a h Mi ) = £. =u k B T\^ (z t - 2^/ 2 z 2 ) + 
O(zf) is the pressure of the an ideal Fermi gas, and b 2 is 
the second virial coefficient which is a function of temper- 
ature only. Using the relations e = Ts + jUfnj + ^x^x — P 
and dP = J2i n id^ + sdT, we have 



A 3 =z m) (l + 2v / 26 2 z i (T)) -2" 
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Z TU)' 



5 P /ZjUj + 
S ~ 2T T 



2V2^f^z, Zl , 



_ + 2V2 —z,z lT — . 



(8) 



(9) 



(10) 



Since b 2 = 1/2, and db 2 /dT = at resonance [9] [3], we 
recover the universal thermodynamics in I. 

It is also useful to use (T, n) instead of (T, fx) as vari- 
ables. The following relations are applicable to all scat- 
tering lengths and can be compared with experiments. 
Eliminating Zi in cq.(8) to (10), we have 

P = k B T(n + [2~ 5/2 (n 2 + n\) - 2V2b 2 n ] n i ]\ z ), (11) 



-k B T 



(n 2 + n?)A 3 
n + ^ ^ n T«l 



(12) 



where (f> = 2 V / 2A 3 (6 2 - lTdb 2 /dT). From eq.(ll) 
and (8), it is easy to calculate isothermal compressibil- 
ity k t ~n^ 1 (dn/dP)r and isothermal spin susceptibil- 
ity Xt~ [dm/dv)T- Their deviations from ideal gas 
values are Ak t = V2b 2 n\ 3 /(nk B T), A Xt = -(n/k B T) 
(5) (\/262«A 3 /4), hence u 2 Ak t /Ax t = -4. It is easy to 
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derive the same results for the Bose gas, which is eq.(ll) 
and (12) with a minus sign in the n 2 terms. 

III. Degenerate Normal Gas: This is the case x = 
k B T/fi « 1. For small spin polarization, the pressure 
can be expanded in x and v/n as P = P(")(T, 



A 3/2 



h^{Bk B Tf 
8(A M )2 



15(Ct/) 2 
32(A/i) 2 



(13) 



where coefficients A, B, C are universal numbers. They 
are written in this form to simplify later discussions. The 
absence of linear T term is due to vanishing entropy at 
T = 0. For ideal Fermi gas, A = B = C = 1. Stability 
conditions ds/dT, dn/dfi > imply that A, B > 0; and 
C > unless the system is ferromagnetic. From eq.(3) 
and the relation to = dPjdv, we have n = n^ n \T, /i, v), 



nAAfj) 



iT 2 (Bk B T) 2 3 fCv 
8(A/i) 2 32 \Afi 



(14) 



and to = [3c/(8A)][n f (Afj,)/(Afi)](Cu), where n f (Afi) = 
A 3 / 2 riy(/i). These two equations for n and to readily 
give the density profile Uj(x) (i =T,|) in non-uniform 
potentials Vi (x) within LDA by substituting pn (x) = Hi — 
Vi(x). Note, however, that these relations are valid only 
for fii(x) » k B T and v « [m. As one approaches the 
surface of the cloud, density decreases and the system 
will switch to Boltzmann regime in surface regions where 
e iH(x)/k B T <<: wu - n densities given by eq.(8). 

To find an accurate formula interpolating between de- 
generate and Boltzmann limit, we note that (consider- 
ing the case v = for simplicity) deep in Boltzmann 
regime, z = e ^/ kBT << 1, eq.(8) is simply n = z/X 3 and 
is the high temperature limit of the ideal gas relation 
n = n id(V' T )i n id(^,T) = X~ 3 f 3/2 (z), where f 3/2 (z) is 
the Fermi integral [10]. On the other hand, eq.(14) is pre- 
cisely the low temperature expansion of the ideal gas rela- 
tion nid(fi,T) with the substitution (p,T) — ► (A/i,BT). 
The desired interpolation will then be of the form 

n (B) (T, ft) = n ld (A(x)n, B(x)T), x = fx/k B T (15) 

where {A{x), B(x)) are functions of x (as required by Uni- 
versality Hypothesis) such that (A(x), B(x)) — > (A, B) as 
x » 1, and (A(x), B(x)) — > (1, 1) as z = e x « 1. Since 
the switching from degenerate to Boltzmann regime take 
place at fi <~ k B T, a simple extrapolation is 



A(x) = 



Ae z 



1 



+ 1 



B(x) = 



Be 1 



1 



+ 1 



(16) 



Eq.(15) and (16) form the desired extrapolation. The 
density profile in a trap calculated within LDA using 
eq.(15) and (16) is shown as the dashed curve in the 
lower figure in figure 1. 



To derive relations related to experiments, we invert 
the relations n = n(T, /i, v), to = m(T, /x, v) to express 
jj,, v, and hence e = 3P/2 in terms of (T,n,m). To the 
lowest order in k B T/£p, we have fi = fi^ n \T, fj,, v), 

,<■» .«•,(«)(! -W)IA, % = (17, 



E W = [3n££(n)/5] [1 + 51U] /A, 



(18) 



where W 



' + ( 



2Am' 
3Cn 



Alternatively, we can 
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write /iW = ££(1+/3J and e W = (3n££/5)(l+/3J. The 
numbers /3 and (3 e are the interaction parameters mea- 
sured in ref. [5] and ref. [6] . We show here that they have 
opposite temperature and spin polarization corrections, 
which differ from each other by a factor of 5. 

IV. Superfluid at resonance: This is the case where 
universal thermodynamics proves very useful. We shall 
consider supcrfluids with zero spin polarization, (hence 
v = 0). Within Ginzburg-Landau theory, the difference 
in grand potential between a superfluid with order pa- 
rameter ^(r) = ( - 0|V'l) an d a normal fluid nrar super- 
fluid transition is fi[tf] - fi(") = / drw[#(r)], 



w[*] = if |W| 2 - r 2 |*| 2 + r 4 |*| 4 /2 



(19) 



where K,r 4 > 0, and r 2 vanishes at transition. The equi- 
librium potential is SI = fi^o], where ^ is the minimum 
of eq.(19). According to Universality Hypothesis, K, r 2 
and r 4 are functions of k B T and /i only. The condition 
for transition r 2 (T = T c ,/i) = implies T c = T c (fi). 
Using dimensional analysis and expanding r 2 , r 4 and K 
about T c , we have k B T c {p) = j(Afi), K = £ > h 2 /{2M); 



and to the lowest order of 1 — T/T c 



a 2/ u(l-T/T c ), 



r 4 = a4/j,/nf(A/j,), where (7, £, a 2 , a 4 ) arc universal 
numbers characterizing the superfluid state near T c . 
When T > T c , the system is normal with $ = 0, and 



P = p(«)( TjjU;0 ), n = nW(r,|i,0). 



(20) 



For T < T c , we have |* D | 2 = r 2 /r 4 , and P = P< n )(T,^) + 
r|/(2r 4 ). Explicitly, we have |* D | 2 = an f (A f i)(l-T/T c ), 
a = a 2 /a 4 , and P = P( s )(T,/i), 

P(^P(«) ( T,,,0) + ^4^(1-^) 2 , (21) 

where D = 5a 2 /(4a 4 ), and x = k B T/fi. Eq.(21) then 
implies n = n^ s ' {T,p), 



n 



s )=n(")(T, Ai ,0)+n / (A M ) J D (l - ^ ^1 



a; 
57 



(22) 



where 7 = 7A, x = k B Tj\i. Eq.(22) and eq.(15) 
provide a simple method to construct the density pro- 
file of a trapped superfluid within LDA: We first plot 



3 



/i(r) = p — V(r) and 7/u(r) as a function of r. (See 
fig.l) The regions where fc^T < jp(r) and 7/i(r) < fc^T 
correspond to superfluid (SF) and normal (N) region. 
The latter is further separated into degenerate normal 
regime (DN), 7/1 >> k B T, and Boltzmann (B) regime, 
e M/feBT << 1 Tng densities m (gp) anc j are gi ven by 

cq.(22) and (15) respectively. The "superfluid bump" in 
fig.l is also obtained in ref. [1]. Here, we show that it is 
a necessary consequence of the Universality Hypothesis 
and display its general structure. 

To express quantities in terms of T and n, we invert 
eq.(22) and then obtain p = p^(T,n), e = e( s )(T,n). 
If superfluid transition takes place in degenerate regime, 
e = (nBk B T c ) 2 /(8(Ap) 2 ) = (nk B -f) 2 /8 << 1. In that 
case, it is simple to show that 

p^ = p^ (T, n, 0) - (£° F (l if) (23) 



e «= e W(T,»,0)-^(l-^) (24) 

where y = k B T/£° F (n), k B T c = jAp = -/Spin), ( = 
2D/(3A) = 5al/(6Aa 4 ). Eq.(24) and (18) imply a uni- 
versal jump [(c s — c n )/c n ] Tr = 12(A/(5it 2 j 2 B 2 ) cross T c . 

In the superfluid phase, a "second" sound mode (112) 
must exist in addition to the first (or ordinary) sound 
U\ . In the col lisional limi t, their velocities ar e u\ — 
y/ (dP/(Mdn) a and u 2 = \J o 2 p s j {p n {do /dT)p) respec- 
tively [12], where a — s/(Mn) is the entropy per unit 
mass, p s is the superfluid mass density which can be 
obtained easily from the gradient term of eq.(19) to 
be Ps = lMn{\ - T/T c ), I = 4£aa [11]; p„ is nor- 
mal fluid mass density, p s + p n = M n. We th en have 
u\ = 2£ F /(3MA), and u 2 / u i = QV 1 - T Q 2 = 

!(75) 2 i{i+ i r(^) 2 }- 1 [is]. 

V. Bosons in unitarity limit: Since Universal- 
ity Hypothesis makes no reference on statistics, eq.(13) 
should also apply to stable Bose system at unitarity, 
{y = for single component gas). However, it provides no 
information about phase coherence and hence does not 
guarantee that the system is a superfluid [14]. Never- 
theless, it predicts that a stable Bose system (superfluid 
or normal) will have a fermionic energy density in the 
unitarity limit, a fact can be tested by experiments. 
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k B T/(<yAp) = 0.65. 



[1] M. Holland, et.al, Phys. Rev. Lett. 87, 120406 (2001); 

M.L. Chiofalo et.al, Phys. Rev. Lett. 88, 090402 (2002) 
[2] H. Heiselberg, Phys. Rev. A. 63, 043606 (2001). 
[3] T.L. Ho and E. Mueller, cond-mat/0306187. 



4 



